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A precise determination of the decay constant of the Dg-meson in 
quenched QCD * 

Andreas Jiittner and Juri Rolf (ALPHA collaboration) ^ 

^Humboldt-Universitat zu Berlin, Institut fiir Physik, Newtonstr. 15, 12489 Berlin 

We present a precise determination of the leptonic decay constant Fd^ of the Dg-meson in quenched lattice 
QCD. This work is particularly focused on the analysis and discussion of all sources of systematic errors. We 
simulate directly at the physical quark masses for five different lattice spacings between 0.1 fm and 0.03 fm using 
0(a)-improvement. The finest lattice is still work in progress. After taking the continuum limit and setting the 
scale with the Kaon decay constant Fk = 160 MeV we arrive at a value of Fo^ ~ 252(9) MeV. Setting the scale 
with the nucleon mass instead leads to a decrease of about 20 MeV of Jb. . 



1. Introduction 

The Dg meson consists of a charm and a strange 
quark. Its direct simulation on the lattice is pos- 
sible without having to use a chiral extrapolation 
for the light quark or an effective theory for the 
heavy quark (e.g. HQET, NRQCD). Therefore 
the Dg is an ideal testing ground as its proper- 
ties are expected to be similar to those of other 
heavy-light systems (e.g. B(s)). 

The decay constant Fd^ parameterizes the lep- 
tonic decay of the Dg meson into a pair of leptons 
with aW 'AS the intermediate vector boson. It is 
defined through the QCD matrix element 

(0|A^(0)|D,(p)) -ip^Fo,, (1) 

where = S7;_j75C is the axial current. 

Fd^ is an important input parameter to CKM- 
physics in phenomenology. Together with a pre- 
cise computation of the decay constant in the 
static approximation jll2j . it can also be used 
to investigate down to which heavy quark mass 
HQET can be applied safely. With this at hand 
one can then eventually get a reliable estimate for 
Fb from an interpolation between the relativistic 
simulation and the static approximation 

Previous quenched results for Fd^ have been 
summarized in 0], which gives a world average 
of 230(14) MeV. Ref. is a recent calculation 
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of Fd^ with an alternative approach which gives 
a value of Fd, = 239(10) MeV. The PDG quotes 
Fd, = 285 ± 19 ± 40 MeV |7j. 

The aim of this project is to obtain a precise 
value for the decay constant Fd^ with a combined 
statistical and systematic error of 3%, apart from 
the quenched approximation. 

The systematic errors can be kept under con- 
trol by using non-perturbative 0(a)-improvement 
j8l9) and renormalization together with a well 
controlled approach to the continuum. In partic- 
ular, we simulated at four different values of the 
lattice spacing between a — 0.1, . . . , 0.05 fm. At 
the moment we simulate with an even finer lat- 
tice (0.03 fm) and we plan to add the results to the 
analysis in the future. We also analyzed the sys- 
tematic error due to the limited solver precision 
and the systematic error due to the contamina- 
tion of the plateaus by excited states and glue- 
balls. Finite volume effects have been shown to 
be negligible at our choice of the lattice size and 
at our choices of the quark masses [TT] . 

2. Strategy and Simulation 

Our simulation parameters are based on previ- 
ous work of the ALPHA-coUaboration. We use 
the Sommer scale ^21 with tq = 0.5 fm to fix 
the lattice spacing. As shown in j^, this choice 
of ro is equivalent to setting the scale with the 
decay constant of the K-meson, Fk — 160 MeV. 
The hopping parameters of the strange quark and 
the charm quark have been taken from publi- 
cations [11113) on the corresponding renormal- 
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Tablel. Simulation data 



p 


^mcas 


L/a 


L/ro 




6.0 


380 


16 


2.98 


4.972(22) 


6.1 


301 


24 


3.79 


4.981(23) 


6.2 


251 


24 


3.26 


5.000(25) 


6.45 


289 


32 


3.06 


5.042(29) 



ization group invariant masses for all values of 
the lattice spacing. In these works, Kg and Kc 
were determined using the masses toTODb = 4.99 
and rgmK — 1.57 of the Dg- and K-meson as 
phenomenological input, while neglecting isospin 
breaking and using Ms/Mught — 24.4 ± 1.5 with 
^flight = ^{Mu + Md) from chiral perturbation 
theory [Tl) . 

All calculations have been performed in 
quenched QCD using clover-improved Wilson- 
Fermions and Schrodinger Functional |15) bound- 
ary conditions. In the time directions all fields 
obey Dirichlet boundary conditions and in the 
space directions they obey periodic boundary 
conditions (up to a phase in the fermion fields, 
which we chose to be 6 = 0.5). The fcrmionic 
boundary fields (,( and C',C' Eire used to con- 
struct the mesonic boundary sources ^B] 

O = ^EUyhsUz), at a;o = 0, 

Starting from the correlation function 

fiixo)^-^(Aiix)O), (3) 

with the improved axial current A^, one derives 
the expression 

FdAxo) ~ -2Za{1 + bAamq){mB,L^y^/^ 
J.I 

X *^ ^(xn-T /2)m-n,. 

X |l - r;2=e-^"^ - jf^e'iT-^o)raG^^ 
+0(a2) 

for the decay constant. A detailed discussion 
can be found in |17| . Here, Za renormalizes 
the axial current and the term linear in 5a sub- 
tracts 0(a)-artefacts proportional to the bare 
subtracted quark mass mg. Both coefficients were 
obtained non-perturbatively on the lattice |10I9| . 



Figure 1. The simulation result for the bare decay 
constant for = 6.0,6.1,6.2,6.45 (from top to 
bottom). The filled symbols mark the a priori 
chosen range of the plateau. 

Plateaus of F„ 




02 , , 1 , , , 1 , , , L 

2 4 6 



The pseudoscalar mass ttid^ is obtained as effec- 
tive mass and (m-Q^L^^^^^'^ normalizes one par- 
ticle states. Ja ^'^'^ /i — — ^(C'C), which can- 
cels out the normalization on the boundary me- 
son sources, have to be determined in the simula- 
tion. As j\ decays exponentially with the pseu- 
doscalar mass, one expects ^^^^^(a^o) to exhibit a 
plateau at intermediate times, where the contri- 
butions r^e~^°^ of the first excited state and the 
contribution 7y<^g-(^-^o)™G from the glue- 
ball both are small. A plateau average can then 
be performed to increase the signal. 

Some simulation parameters are summarized in 
table 1. We set T = 2L. All simulations were 
done at approximately constant physics, which is 
important for the scaling study. Fig. 1 shows the 
simulation result for the unrenoramlized Fd^{xo). 

3. Error analysis and results 

The simulations for the first four values of /3 
were done in single precision arithmetics on the 
APE-machine at DESY Zeuthen. As we are deal- 
ing with heavy quark propagators, we have to 
check for rounding errors which may arise during 
the calculation. We also want to guarantee suffi- 
cient solver convergence for the propagators. The 
corresponding checks were done with an adapted 
version of the MILC code ^ on the IBM p690 
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Figure 2. Continuum extrapolation for Fd^ 




computers at the HLRN ,5_. Starting from the 
same field configuration, we calculated the prop- 
agators in single and double precision arithmetics 
and with different solver residuals, ft turned out 
that the impact of these changes on Fd^ is below 
1 per mil. 

For the definition of plateau ranges of |0J 
where the relative contribution of excited states 
and glueballs is below a threshold of 5 per 
mil, one needs estimates of A and ma- These 
were obtained self-consistently from linear fits to 
log(^barc(^^)_^piatoau^^ whcrc fP'^*""" is the av- 
erage of an a priori chosen plateau (cf. fig. 
and for details). We found that a suitable 
range for the plateau is xq = 4ro, . . . , 5ro for all 
values of (3. 

Fig. 10) shows the corresponding plateau aver- 
ages for all simulated values of (3. We excluded 
the coarsest lattice in the continuum extrapola- 
tion which was done linear in (a/ro)'^ because 
of 0(a)-improvement. With &a taken from [S] 
we finally quote roFo^ = 0.638(24) as the main 
result. We also used 1-loop perturbation the- 
ory for 6a since taking 6a from |H] involves 
an extrapolation of the data. One then gets 
''o-Fb, = 0.631(24). Excluding the two coarsest 
lattices instead leads to roFu^ — 0.630(34). Fi- 
nally, using ro — 0.5 fm, our main result in phys- 
ical units is Fd, = 252(9) MeV. When using the 
nucleon mass to set the scale, which corresponds 
to taking tq = 0.55 fm, we found that .Fba de- 
creases by 20 MeV 01!. 

4. Conclusion 

We present a direct calculation of the heavy- 
light decay constant Fd^ with the final result 
Fd, = 252(9) MeV. The error matches that of 
future experiments, e.g. CLEO-c. We aim at ex- 



tending the scaling study in the future and these 
calculations are under way. We supplement this 
analysis with more data around the charm mass 
ini ■ We estimate that Fb^ decreases by by 20 MeV 
under a scale shift of 10%. 
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